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ABSTRACT 

We have investigated the consequences of merging double white dwarf systems by 
calculating evolutionary models of accreting white dwarfs. We have considered two cases; 
a massive C-0 white dwarf of ~ 1M Q accreting C-0 mixture, and a low mass white dwarf 
with an initial mass of 0.4M Q accreting matter composed mostly of helium. 

The accretion rate of the C-0 white dwarf is assumed to be 1 x lO _5 M y~ 1 . After 
carbon burning is ignited at M r ~ 1.04M Q , the flame propagates inward due to heat 
conduction. By inserting enough grid points to resolve the structure of the flame, we have 
obtained almost steady burning in most phase of evolution, but we have found a new 
phenomenon that the strength of flame sometimes oscillates due to a thermal instability 
in early phase of the evolution. In calculating evolutionary models, we have occasionally 
employed a steady-burning approximation, in which the propagation speed of flame is 
given a priori. We have considered two extreme cases for the interior abundance of the 
massive white dwarf: Xc = 0.5 and Xc = 0.2. For Xc = 0.5, the flame becomes very 
weak at M r ~ O.4M and the inward propagation stall there. But a few thousand years 
later, the flame is reactivated due to contraction and propagates to the center. For 
Xc = 0.2, the flame propagates smoothly and reaches the center in ~ 1000 years. In 
both cases the C-0 mixture has been burned into an O-Ne-Mg mixture without causing 
an explosive phenomenon. 

For helium accreting low mass white dwarfs, we have considered accretion rates of 
1 x 10~ 7 and 1 x lO _6 M y _1 . After a fraction of M Q is accreted to the white dwarf, 
helium is ignited in the outer part and a shell flash occurs. Such a shell flash diminishes 
when about 10% of helium in the convective shell is burned into carbon and oxygen. 
The next shell flash occurs at a shell interior to the previously flashed shell. After 
less than 30 shell flashes, the helium ignition occurs at the center and steady burning 
begins. Thus, the merging produces a helium star which burns helium at the center. 
The mechanism of the propagation of the burning shell in this case is compressional 



heating during inter-pulse phases, in contrast to the case of carbon-burning flame where 
the conduction drives the inward propagation of the flame. We infer that such a low 
mass double white dwarf system could be a progenitor of the AM CVn stars. 

Subject headings: stars: evolution - stars: white dwarfs - stars: binary - stars: chemi- 
cally peculiar - stars: supernova 
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1. Introduction 

It is expected theoretically that a fraction of close 
binary systems become double white dwarf systems 
(Iben & Tutukov 1984; Webbink 1984). Obscrvation- 
ary, several double white dwarf systems have been 
found by Nather, Robinson, & Stover (1981), Saffer, 
Liebert, & Olszewski (1988), Bragaglia, et al. (1990), 
Marsh, Dhillon, & Duck (1995), and Marsh (1995). 
Although low mass white dwarfs (M < 0.45M Q ) are 
often found in binary systems (Marsh et al. 1995), 
no massive double white dwarf system with a short 
period (less than several hours) has been found de- 
spite several searches by e.g., Robinson & Shatter 
(1987), Foss, Wade & Green (1991), and Bragaglia 
et al. (1990). Bragaglia (1995) discussed the present 
status of the search and possible explanations for the 
null result. In this paper, we assume that there exist 
some massive double white dwarf systems as well as 
low mass systems. 

Double white dwarf systems having short enough 
periods are expected to merge within the Hubble 
time by the effect of gravitational wave radiation. 
Since the radius of a white dwarf is proportional to 
M -1 / 3 , once the less massive component fills its crit- 
ical Roche lobe, a dynamical mass transfer occurs if 
the secondary to primary mass ratio is sufficiently 
large (Cameron & Iben 1986; Benz et al. 1990). 
The consequence of the rapid mass transfer seems to 
be the formation of a hot thick disk around the pri- 
mary white dwarf (e.g., Hachisu, Eriguchi, & Nomoto 
1986; Mochkovich & Livio 1990). At least a part of 
the disk matter eventually falls on the white dwarf to 
complete the merging process. Iben & Tutukov (1984) 
and Webbink (1984) argued that if the total mass of 
the two white dwarfs is larger than Chandrasekhar's 
limiting mass, the merging would lead to a Type la su- 
pernova. Whether it occurs or not, however, depends 
on the accretion rate onto the white dwarf, which, 
in turn, depends on the timescale on which the an- 
gular momentum distribution changes significantly in 
the thick disk. As discussed by Mochkovich & Livio 
(1990), the timescale can be as short as a few days if 
the thick disk is turbulent. In such a case, the accre- 
tion onto the white dwarf is limited by the Eddington 
limit (M ~ 2 x lO^Moy" 1 for M ~ 1M ) due to the 
effect of radiation pressure. For such a high accretion 
rate, Nomoto & Iben (1985) have shown that carbon- 
burning is ignited in an outer shell of the white dwarf 
due to the effect of comprcssional heating. Saio & 



Nomoto (1985) and Woosley & Weaver (1986) have 
shown that carbon-burning ignited near the surface 
propagates to the center to burn most of carbon with- 
out causing an explosive phenomenon. In these calcu- 
lations the propagation of carbon-burning was simu- 
lated as successive phases of a shell-flash and an inter 
shell-flash, because the burning region (flame) was not 
spatially resolved. 

To overcome the poor resolution in stellar evolution 
models, Timmes, Woosley, & Taam (1994) calculated 
the propagation of flame in initially homogeneous 
medium with fixed boundary temperatures. They re- 
solved the structure of flame by using an adaptive 
mesh point method. They obtained the flame speed as 
a function of the hot side boundary temperature and 
the cooler side density. Combining the result with the 
assumption of balanced power between the integrated 
nuclear energy generation rate and neutrino emission 
rate, they have concluded that the flame will reach 
the center if the mass of the white dwarf is greater 
than ~ 0.8M Q . On the other hand, Kawai, Saio, & 
Nomoto (1987) have shown that carbon burning can 
ignite in the outer part of an rapidly mass-accreting 
white dwarf if it is more massive than l.O7M . These 
results indicate that once carbon-burning is ignited 
the flame propagates throughout the interior of the 
white dwarf to convert the C-0 mixture into an O- 
Ne-Mg mixture. Therefore, if the merging of a double 
white dwarf system leads to an accretion rapid enough 
to cause an off-center carbon ignition in the primary 
white dwarf (M £ 2.7 x lO _6 M y _1 ; Kawai et al. 
1987), no Type la supernova is expected. Since it 
is very important to know what is the final product 
of the merging of double white dwarfs, we have re- 
investigated the flame propagation in a white dwarf 
by inserting enough grid points to resolve the flame 
structure. 

In a helium white dwarf accreting helium at a rate 
faster than 2 x lO _8 M0y _1 helium burning is ignited 
off-center (Nomoto & Sugimoto 1977). When we 
made computations for the massive C-0 white dwarf 
before, we also calculated helium white dwarf mod- 
els. We found that helium burning shell propagates 
inward to the center by repeating mild flashes, and 
that the white dwarf is converted to a helium-main- 
scquence star. Some of the results have been pre- 
sented in Nomoto & Hashimoto (1986, 1987). Iben 
(1990) investigated similar cases and claimed that the 
products of merging low mass white dwarfs are sdB 
and sdO stars (sec, e.g., Hebcr (1991) for a review 
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of sdB and sdO stars). In our present paper we also 
presents results of additional computations for helium 
accreting low mass helium white dwarfs. 

2. C-O white dwarf 
2.1. Structure of flame 

In our previous calculation (Saio & Nomoto 1985), 
the inward propagation of the carbon-burning shell 
appeared as a series of successive shell flashes. This 
is an artifact brought by insufficient spatial resolu- 
tion. The actual evolution should be a nearly steady 
progression of the burning front or flame (Timmes et 
al. 1994). In this paper, we present the result of 
our re-computations for this problem done by insert- 
ing enough grid points to resolve the carbon burning 
flame. Input physics and the initial model are the 
same as m Saio k Nomoto (1985). Figure shows 
the flame structure at selected phases when the flame 
is located at M r ~ O.6M . The flame, whose center is 
defined by the peak of nuclear energy generation rate 
per unit mass e n , moves inward keeping the neigh- 
boring structure similar. Just interior to the flame 
there is a thermal front having an extremely steep 
temperature gradient. The front moves inward due 
to the conductive energy flow. Matter experiences a 
large heating by the passage of the front. The density 
decreases appreciably between the thermal front and 
the flame center. 

Since the carbon abundance decreases outward in 
the flame, the maximum temperature occurs at a zone 
exterior to the flame center. The position of the max- 
imum temperature is recognized in Figure |l| as the 
point where the direction of energy flow changes (i.e., 
the sign of L r changes; dotted lines are used if L r < 
and long dashed lines if L r > in Fig. |l|). Figure 
[l] indicates that most of the energy produced in the 
flame flows inward and is used to move the thermal 
front inward. 

Behind the flame (i.e., exterior to the flame) a shell 
convection zone develops. The bottom of the con- 
vection zone is recognized in Figure |l| as sudden de- 
creases of Xq and e n . The convection zone at these 
stages extends to M r ~ 0.73A/ Q . A considerable 
amount of energy is generated from carbon-burning 
in the convection zone. For the model in the mid- 
dle panel of Fig. |, L n = 7.64 x 10 6 i o , in which 
4.72 x 10 6 L Q is produced in the flame, and the re- 
maining part is produced in the convective zone, and 
L y = 3.06 x 10 6 Lq. Around the outer boundary of 



the convection zone the outward energy flow is re- 
duced to L r ~ 1.8 x 10 4 -Lq due to the neutrino emis- 
sion in the convection zone. Roughly speaking, 40% 
of the energy generated by carbon burning is lost by 
neutrino emission and the remaining energy is con- 
ducted inward to heat matter between the flame and 
the thermal front. Nomoto (1984,1987) and Garcia- 
Bcrro, Ritossa, & Iben (1997) have obtained similar 
results for a carbon-burning flame in the asymptotic 
giant branch phase of a 9M Q model. 

2.2. Propagation of flame 

Keeping grid spacing fine enough to resolve the 
thermal front and the flame demands very small time- 
steps to follow the flame propagation. To alleviate 
the computational effort we used a steady burning 
approximation occasionally in evolutionary computa- 
tions. In the steady burning approximation, we do 
not resolve the structure around the flame; i.e., the 
thermal front and the flame is represented by a grid 
point. The total nuclear energy generation rate L n is 
calculated by 

L n = X c Q\dm { /dt\, (1) 

where Q is the energy yield from carbon burning, for 
which we have adopted Q = 4.4 x 10 17 ergs s _1 g _1 , 
and mf is the mass coordinate of the flame center. 
The local energy generation rate per unit mass e n 
was set to be zero except at the grid point corre- 
sponding to the flame; i.e., e n j = if j ^ J and 
e n j — L n /{M r j + i — mf), where Jth grid point cor- 
responds to the flame and M r j represents the mass 
coordinate at the jth grid point. The value of rhf 
for the steady burning calculations was obtained from 
a few thousand preceding models obtained by usual 
evolution calculations with fine grid spacings. 

In our present investigation, evolutionary calcula- 
tions were started with a model from Saio & Nomoto 
(1985). The model has a carbon-burning flame at 
mf ~ 0.85Mq. Interior to the flame, carbon and oxy- 
gen abundances are (Ac,Ao)=(0.5,0.5). To see the 
effect of the carbon abundance in the interior, we have 
also computed a set of evolutionary models starting 
with the initial model in which the abundance inte- 
rior to the flame had been artificially converted to 
(Ac,Ao)=(0.2,0.8). In both cases, the accretion rate 
was kept constant at 10 _5 M Q y _1 . 

Figures || and ^ show evolutionary changes in the 
nuclear energy generation rate L n (solid line) and the 
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neutrino energy loss rate L u (broken line) as a func- 
tion of the mass coordinate of the carbon-burning 
flame mf for Xq — 0.5 and 0.2, respectively. In these 
figures, the phases with constant L n correspond to the 
evolution phases during which the steady burning ap- 
proximation was applied. As can be seen from these 
figures, one set of calculations with the steady burn- 
ing approximation was stopped when mf changed by 
~ 0.05M Q . Large fluctuations in L n at the transi- 
tions from the steady burning calculation to the nor- 
mal evolutionary calculation are artificial phenomena 
caused by inserting grid points just interior to the 
flame. 

In most phase of the progression of the carbon- 
burning flame the nuclear energy generation rate L n 
exceeds the neutrino energy loss rate L v . As discussed 
above, the difference is approximately the amount 
of energy conducted inward. This energy balance is 
similar to that found by Garcia-Berro et al. (1997) 
but somewhat different from the balanced power as- 
sumption employed by Timmes et al. (1994). For 
Xq = 0.5, the inward propagation of the flame stalls 
at M r ~ 0.4M Q , which appears as a sharp dip of 
L n in Fig. ||. We discuss this phenomenon more 
below. In both cases the flame reaches to the cen- 
ter and exhausts carbon interior to the first ignition 
shell, which confirms our previous conclusion (Saio & 
Nomoto 1985). 

Figures |] and || show the time variation of the dis- 
tributions of temperature and of density, respectively, 
during the flame propagation for Xq = 0.5. The flame 
temperature is maximum when the burning front is 
located at M r ~ 0.8M©, i.e., when m f ~ 0.8M Q . 
Then it decreases gradually as the flame propagates 
toward the center. Since pressure is continuous across 
the flame, density decreases as matter gets into the 
flame. The density just interior to the flame is max- 
imum when mf ~ 0.9M Q . It decreases as the flame 
propagates inward until it becomes minimum when 
mf ~ O.5M , and then increases slowly as the flame 
approaches the center. These variations depend on 
the flame speed, bulk expansion rate of the white 
dwarf, density gradient, the degree of electron degen- 
eracy, etc. 

The central temperature and density decrease as 
the flame propagates. This variation is caused by the 
net gain of energy because the nuclear energy genera- 
tion rate exceeds the neutrino loss rate in most phases 
of the flame propagation. (When the neutrino loss 
rate exceeds the nuclear energy generation rate, both 



the central temperature and density increase [see be- 
low].) The decrease of the central temperature seen 
in Fig. [I] is larger than the case shown in Saio & 
Nomoto (1985). This can be explained as follows: In 
Saio & Nomoto (1985), the spatial resolution was not 
enough so that the flame propagation appeared as a 
series of shell-flash and inter shell-flash phases. Dur- 
ing an inter shell-flash phase, the neutrino loss rate 
exceeded the nuclear energy generation rate. There- 
fore, the net gain of energy in the present calculation 
is higher than the case in Saio & Nomoto (1985). 

The time variation of m,f for Xq = 0.5 is shown in 
Figure ^ by the dashed line. The solid line and the 
dotted line in this figure show the position of ther- 
mal front for Xq — 0.5 and 0.2, respectively. The 
shell having logT = 7.7 well represents the position 
of the thermal front, which is located just interior to 
the flame (see Fig. [l]). The propagation speed of the 
flame for Ac = 0.5 decreases significantly and the 
flame stalls for a while at M r ~ O.4M . This phase 
corresponds to a deep dip in L n seen in Fig. ^[ Even 
when the flame propagation stops, the thermal front 
gradually moves inward but the steepness of the tem- 
perature gradient decreases (Fig. [j] below). After the 
flame have stalled there for ~ 4000 years, it becomes 
active again and propagates toward the center at a 
speed considerably slower than before. After the re- 
activation, the carbon abundance in the convective 
shell is ~ 0.1, which is considerably higher than the 
case before the flame stalled. After the flame reached 
the center, the remaining carbon has been burnt qui- 
etly in a few thousand years. 

The flame for Xq — 0.2, on the other hand, reaches 
the center smoothly in ~ 1000 years and the thermal 
front keeps very close to the flame during the prop- 
agation. Since the flame remains strong throughout 
the propagation in this case, the carbon abundance 
in the convective zone never exceed 0.005. Therefore, 
when the flame reached the center, most of carbon in 
the region interior to the first ignition shell has been 
exhausted. 

Figures [7] and |8| show variations of temperature and 
density distributions, respectively, in the phase dur- 
ing which the propagation of the flame almost stopped 
and was re-activated later for Xq = 0.5. Upper and 
lower panels of these figures show evolution before and 
after the re-activation of flame, respectively. The di- 
rection of evolution is recognized by the inward prop- 
agation of the thermal front. The last lines in the 
upper panels are the same phase as the first lines in 
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the lower panels. Before the re-activation, as the ther- 
mal front propagates inward, the temperature gradi- 
ent at the front becomes less steep and the density 
jump gets smaller and less apparent. During this 
phase, the flame becomes weak and hardly moves at 
M r ~ 0.42M© (Fig. |). Since L D < L v in the phases 
shown in the upper panels of these figures, the density 
interior to the front increases with time. The burn- 
ing becomes active again at M r ~ 0.4M Q , when the 
thermal front is located around M r ~ O.35M . A 
new thermal front backed by the activated flame ap- 
pears outside of the old front and propagates inward 
to engulf the old one at M r ~ 0.345M Q . 

The time spent for the flame to reach the center is 
shorter for Xq = 0.2 than the case of Xq = 0.5. 
Moreover, for Xq = 0.5 the flame in the present 
calculation takes longer to reach the center than in 
our previous calculation (Saio & Nomoto 1985), in 
which it took about 5000 years for the flame to prop- 
agate through the interior. Such differences can be ex- 
plained by the fact that the net-energy gain (L n — L v ) 
in the present Xq = 0.5 case is larger than the 
Xc = 0.2 case or our previous calculation. A larger 
energy-gain makes gentle the temperature increase in 
the inner part of the flame, because a larger energy- 
gain decreases more the density at the ignition tem- 
perature, thus making the degree of degeneracy there 
lower. A lower temperature gradient makes the prop- 
agation speed lower. 

2.3. Discussion 

Timmes et al. (1994) evaluated propagation speeds 
of carbon-burning flames in a isobaric region which 
has initially homogeneous temperature and density. 
They obtained flame speed in such systems as a func- 
tion of the boundary temperature of the hotter side 
and the density of the cooler side (initial density) of 
the flame. The boundary temperature in the cooler 
side (i.e., initial temperature) was fixed at 5 x 10 8 K, 
which is much higher than in our models. Figure ^ 
shows the propagation speed \dmt/dt\ of the flame de- 
rived from the actual propagation in our evolutionary 
models for Xc = 0.2 (filled circles) and that taken 
from the table of propagation speed for Xq — 0.2 in 
Timmes et al. (1994) for the density and tempera- 
ture at the flame in our models (open circles). This 
figure indicates that the propagation speed based on 
Timmes et al.'s result tends to be higher than that 
in our evolutionary models. One reason for this dif- 
ference may be ascribed to the fact that the cooler 



side temperature of our models is much lower than 
that in Timmes et al.'s analyses. We note here that 
evolution calculations using the steady state approx- 
imation with the propagation speed of Timmes et al. 
failed before the flame reaches the center because the 
density became too low. The rapid density decrease 
occurred because the energy production rate was too 
high, which indicates that their propagation speed 
seems too high for our white dwarf models. 

We have found a new phenomenon that the prop- 
agation of the flame tends to be oscillatory when the 
peak temperature is relatively high. An example is 
shown in Figs. [lO] and O. The lower panel of Fig. 10 
shows runs of temperature in the interior for every 20 
evolutionary models, while the upper panel shows nu- 
clear energy generation rate L n and neutrino emission 
rate L u as functions of rri{ (time advances from right 
to left of the diagram) . Figure |ll| sh ows L n and L v as 
functions of time. As seen in Fig |l0|, the temperature 
gradient between the ignition point (logT ~ 8.8) and 
the position of maximum temperature becomes gen- 
tle occasionally, and the energy generation rate dips. 
The time variation of L n is characterized by semi- 
regular oscillations with sharp peaks and relatively 
gentle minimums. The typical 'period' is an order of 
O.ly. In contrast to the variation of L n , L v hardly 
changes. 

Such oscillatory behavior of the flame tends to oc- 
cur when the thermal front is resolved by finer grid 
spacings. Since finer grid spacings make the max- 
imum temperature steepness of the thermal front 
larger, we infer that a steep temperature gradient 
plays an important role in the instability. We may 
explain the cause of the fluctuation as follows: When 
the temperature gradient becomes too steep, the in- 
ward energy flow exceeds the energy generation by 
nuclear burning to decrease the flame temperature 
and hence L n . This decreases the steepness of tem- 
perature distribution, which in turn slows the inward 
energy flow. Then, the flame temperature and L n 
increase to recover the previous levels. 

It is known that a combustion flame front in a 
solid medium (in which the diffusion rate of the fuel 
is much slower than the rate of thermal conduction) 
is thermally unstable if the temperature sensitivity 
of the energy generation rate is sufficiently high (see 
e.g., Zeldovich et al 1985). The high temperature 
sensitivity of energy generation rate makes the flame 
thin with a steep temperature gradient. We believe 
that the instability is the same as that of the carbon- 
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burning flame in our accreting white dwarf because 
both of them occur under the condition of a steep 
temperature gradient and a high temperature sensi- 
tivity of energy generation rate. 

3. Heium white dwarf 

The starting model adopted for a helium-accreting 
white dwarf is a chemically homogeneous (Y=0.98, 
Z=0.02) white dwarf of 0.4M©. The white dwarf is 
assumed to be the massive component of a double 
white dwarf system. When the less massive compo- 
nent fills the critical Roche lobe (because of the grav- 
itational wave radiation), the mass transfer is unsta- 
ble if the mass ratio exceeds 2/3. In this case the 
massive component is expected to accrete matter at 
a rate of the Eddington limit, which is of the or- 
der of lO _5 M0yr~ 1 . If the mass ratio is less than 
2/3, stable mass transfer occurs at a rate of e.g., 
3 x KrV" 1 for Mi = 0.2M Q , and 5 x lO^A/eyr" 1 
for M% = 0.1M©, where Mi indicates the mass of the 
less massive component. Iben (1990) investigated the 
consequences of rapid merger processes in which mass 
is accreted at a rate faster than the critical Eddington 
limit. We investigate here cases with mass accretion 
rates lO _7 M yr _1 and 10 -6 MQyr -1 , which are much 
smaller than the critical Eddington limit. 

The chemical composition of accreting matter is 
set to be (Y, Z) = (0.98,0.02) assuming that the less 
massive component is also a helium white dwarf. 

3.1. Results 

Accretion causes heating in the envelope of the 
white dwarf. When a sufficient mass is accreted, the 
triple alpha reaction is ignited. For higher accretion 
rate the total mass at the helium ignition is lower 
and the mass between the ignition shell and the sur- 
face is smaller. The ignition occurs at M r ~ O.44M 
when the total mass becomes M ~ O.52M for M = 
10~ 7 M y _1 ; these quantities are, respectively, M r ~ 
0.41 M Q and M ~ 0.45M© for M = lO^A^y -1 . 
The ignition of the triple alpha reaction leads to a 
shell flash (thermal pulse). 

After the shell flash is over, another shell flash oc- 
curs at a slightly inner shell than the previous one. 
Such shell flashes repeat ~ 20 to 30 times until the 
ignition shell reaches the center of the star. The num- 
ber of shell flashes is larger than for the cases inves- 
tigated by Iben (1990). The difference is attributed 
to the fact that the first ignition shell in Iben's mod- 



els (M r ~ 0.25Mq) is closer to the center than our 
models. Figure [l2| shows the time variation of the 
integrated nuclear energy generation rate L n in so- 
lar units for M = lO -6 M y -1 . The strength of a 
shell flash decreases as the ignition shell moves in- 
ward. The evolution of repeating mild flashes is simi- 
lar to the evolution which occurs after the ignition of 
helium core flash in a red giant star (e.g., Mengel & 
Sweigart 1981). 

The first shell flash is stronger for the lower ac- 
cretion rate, because the density at the ignition shell 
is higher. At the peak of the first shell flash logL n 
reaches up to ~ 7.5 for M = lO _6 M y _1 , while 

- 10.3 for M = lO^Moy -1 . The strength of the 
second and the later shell flashes hardly depends on 
the accretion rate. 

Figure shows the time variation of the mass 
coordinate rri{ where the nuclear energy generation 
rate per unit mass is maximum. The zero point of 
time in this figure is set at the first ignition of he- 
lium burning. The helium burning shell reaches the 
center in ~ 5.7 x 10 5 y for M = lO _6 M y -1 and in 

- 1.2 x 10 6 y for M = lO^A^y" 1 after ignition. The 
inward move of mi is faster for higher accretion rate 
because the total mass is larger. The total mass in- 
creases to about IMq before the burning shell reaches 
the center for M = lu- 6 M y _1 , while O.63M for 
M = lO _7 Af y _1 . At each shell flashes, a convec- 
tion zone appears above the burning shell. One shell 
flash ends when the helium abundance in the convec- 
tion zone decreases to ~ 0.90. Therefore, when the 
burning shell reaches the center, the star becomes to 
have a structure similar to that of the zero-age helium 
main-sequence. 

The inward propagation of the burning shell does 
not stop even if the accretion stops. Figure [l3| in- 
cludes two cases in which accretion has stopped be- 
fore the burning shell reaches the center. In one case, 
A6S45, the accretion (M = 10" 6 Af Q y~ 1 ) has stopped 
when the total mass became O.45M (just after the 
first shell flash) and in another case, A6S50, the accre- 
tion has stopped when the total mass became O.5M 
(when the burning shell was located at M r ~ O.3M ). 
The time in which helium burning reaches the center 
for A6S45 is about 2.7 x 10 6 , which is longer than for 
A6S50 and is similar to Iben's (1990) O.46M case. 
Since the accretion process assumed in Iben (1990) is 
quite different from ours, we can conclude that the 
time is mainly determined by the total mass and is 
insensitive to the accretion rate itself. 
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Figure |l4] shows the evolutionary paths of our ac- 
creting helium white dwarfs in the HR diagram. Af- 
ter accretion starts the luminosity increases along the 
white dwarf sequence until the triple-alpha reaction 
is ignited. Once the ignition occurs, the luminos- 
ity decreases at first and increases later. By the 
end of the first shell flash, the star moves into the 
super-giant region. From this behavior we can in- 
fer that even if the mass transfer is stable and M is 
as low as 1 x 10 _7 MQy _1 , once helium burning is 
ignited, the radius increases to engulf the compan- 
ion. As the burning shell moves inward with repeat- 
ing shell flashes, the star gradually evolves blueward 
to the hclium-main-sequence. The redward extension 
of the evolutionary paths during the phase of repeat- 
ing flashes is larger than Iben's (f990) models, be- 
cause in the latter models the mass above the first 
ignition shell (> O.f Mq) is larger than in our models 
(~ 0.04M Q for M = 1 x 1CT 6 , and - 0.07M Q for 
M = 1 x fCT 7 ). 

When the burning shell reached the center, the 
accretion rate was set to be zero. Then, the total 
mass is 0.63M Q for M = 1 x 1CT 7 and 1.01M Q for 
M = 1 x 10~ 6 . After the burning shell has reached 
the center, the star spends core helium burning phase 
as a helium-main-sequence star which lasts of or- 
der of 10 7 y. The star evolves off the helium main- 
sequence when the helium is exhausted in the core. 
The star with 0.63M Q (dotted line) evolves to be- 
come brighter and bluer and will evolve along with the 
white dwarf cooling sequence. The star with f.0fM Q 
evolves to super-giant region again; the behavior is 
known from the helium star evolution calculations by 
e.g., Paczyhski (1971) and Saio (1995). 

3.2. Discussion 

3.2.1. Comparison with the case of carbon-oxygen 
white dwarfs 

Figure [l5] shows runs of the temperature and the 
rate of gravitational energy release e g around the 
burning shell at selected stages from the peak of a 
shell flash to the next one. The number for each 
curve indicates the corresponding model number. At 
the peak of a shell flash a sharp peak appears in the 
temperature distribution (#95f and #1251), and a 
zone exterior to the temperature peak has negative 
e g , which indicates that it is expanded by absorbing 
heat released by the helium burning. This figure also 
shows that during the inter-flash period, the layer be- 



neath the flashed shell is heated by compression (i.e., 
e g > 0). The next shell flash occurs at a shell in 
the heated layer. Thus, the main driving force of the 
inward propagation of the burning shell is the com- 
pressional heating during the inter-flash phase. This 
differs from the case of the carbon burning flame. In 
the latter discussed in the previous section, 

inward conductive heat flow is the main cause for the 
flame propagation. 

In the case of carbon burning flame the burning 
shell exhausts carbon during the inward propagation 
and most of the released nuclear energy is conducted 
inward or lost by neutrino emission in the convection 
zone just above the flame. In the case of helium burn- 
ing shell, on the other hand, only a small fraction of 
helium is converted to carbon and oxygen during the 
inward progression of the burning shell and neutrino 
emission rate is always much smaller than the nuclear 
energy release. 

3.2.2. Relation with observed double white dwarf 
systems 

Several double white dwarf systems which consist 
of low mass (<; 0.4M Q ) white dwarfs have been found 
(Saffer, Liebert, & Olszewski 1988, Bragaglia, Greg- 
gio, & Renzini 1990, Marsh, Dillon, & Duck 1995, 
Marsh 1995). The secondary components of these 
systems seem to have comparable masses. Since the 
masses of those white dwarfs are too small to ignite 
the triple alpha reaction by the effect of their own 
gravity, these are probably helium white dwarfs. Al- 
though all the systems are detached, some of the sys- 
tems have orbital periods of a few hours and are ex- 
pected to become semi-detached systems in a few bil- 
lion years (Marsh et al. 1995, Marsh 1995). 

When the less massive component of a double 
white dwarf system exceeds its critical Roche lobe, 
mass transfer starts. If the mass ratio of the sec- 
ondary to the primary is higher than about 2/3, mass 
transfer on the dynamical time scale is expected to 
occur (at least in the beginning) to form a common 
envelope, and the primary accretes matter at a rate 
of the critical Eddington limit. Even if the mass ra- 
tio is smaller and mass transfer is stable, a common 
envelope is expected to be formed because the radius 
of the primary becomes the supergiant size after the 
triple alpha reaction is ignited in the primary as dis- 
cussed above (Fig. [l4|) . 

We do not know exactly what is the consequence 
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of the common envelope phase. We infer that matter 
is peeled off from the secondary. If the system loses 
angular momentum effectively enough, the separation 
between the primary and the secondary decreases to 
merge (i.e., 'spiral- in' occurs). The system would be- 
come a single helium star with a mass of ^ O.8M , 
which will eventually become a DB white dwarf. 

If the systemic loss of angular momentum is not 
large, the separation can increase as the mass trans- 
fer proceeds. In this case, common envelope phase 
ends when the radius of the primary shrinks and the 
remnant of the secondary gets out of the primary en- 
velope. Considering the fact that there exist recurrent 
novae, in which in each nova explosion the secondary 
component is engulfed in the expanding envelope, we 
infer that the secondary in the double white dwarf 
system can survive the common envelope phase. In 
the double white dwarf system, however, the sepa- 
ration is much smaller and the duration of the com- 
mon envelope phase is much longer than in the case 
of novae. Therefore, it is expected that a significant 
amount of mass is peeled off from the secondary. If 
the secondary survives, the remaining system would 
have a primary in the helium main sequence phase 
with a mass less than 0.8M Q and a secondary white 
dwarf having a very small mass. After helium is ex- 
hausted in the primary, the system becomes a double 
white dwarf system with a small mass ratio. Such 
a system loses its angular momentum by the gravi- 
tational wave emission and eventually the less mas- 
sive component fills its critical Roche lobe to resume 
mass transfer. (Because of the small mass ratio, the 
mass transfer is dynamically stable.) Such interact- 
ing double white dwarf systems can be identified as 
AM CVn type stars, which are hydrogen-deficient cat- 
aclysmic variables composed of double white dwarfs 
with a small mass ratio (M2/M1 < 0.1; see Warner 
(1995) for a thorough review on the AM CVn stars). 

4. Conclusion 

We have investigated the consequences of merg- 
ing double white dwarf systems by calculating evolu- 
tionary models for accreting white dwarfs. We have 
considered two cases; a massive C-0 white dwarf of 
~ M© accreting C-0 mixture, and a low mass white 
dwarf with a initial mass of OAMq. 

For the C-0 white dwarf an accretion rate of 
1 x 10 -5 MQy -1 was assumed. After the carbon burn- 
ing is ignited at M r ~ l.O4M , the flame propagates 



inward due to heat conduction. We have considered 
two extreme cases for the interior abundance of the 
massive white dwarf: Xq =0.5 and Xq =0.2. For 
Xq = 0.5, the flame becomes very weak at M r ~ 
0.4M Q and the inward propagation stalls there. But 
a few thousand years later, the flame is reactivated 
due to contraction and propagates to the center. For 
Xq = 0.2, the flame propagates smoothly and reaches 
the center in ~ 1000 years. In both cases, most of the 
interior of the white dwarf has been converted from 
the C-0 mixture into an O-Ne-Mg mixture without 
causing an explosive phenomenon. Therefore, if the 
merging of a double C-0 white dwarf system leads 
to an accretion onto the primary white dwarf rapid 
enough to ignite carbon off-center, we cannot expect a 
Type la supernova but the formation of a neutron star 
(Saio & Nomoto 1985; Nomoto & Kondo 1991). We 
have found a new phenomenon that the flame propa- 
gation tends to be oscillatory when the peak temper- 
ature is relatively high. 

For the helium white dwarf we have considered ac- 
cretion rates 10 _7 M Q y _1 and lO _6 M0y _1 . In each 
case, when a certain amount of helium is accreted, 
the triple-alpha reaction is ignited in the outer part 
and a shell flash (thermal pulse) occurs. The helium 
burning shell moves inward by repeating alternatively 
shell flashes and inter-flash compression. The burn- 
ing shell reaches the center on a timescale of order 
10 6 y after the first ignition. The inward progression 
docs not stop even if the accretion has stopped after 
the first shell flash. At each shell flash a convection 
zone is formed. When the helium abundance in the 
convective layer decreases by ~ 10%, the shell flash 
diminishes. When the burning shell reached the cen- 
ter the star has a structure similar to that of zero-age 
helium main-sequence star. Such a star is expected 
to evolve as a helium star to finally become a DB 
white dwarf. We have also discussed the possibility 
that a low mass double white dwarf system could be 
a progenitor of the AM CVn stars. 
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Fig. 1. — Runs of density p, temperature T, carbon 2 °" m 06 8 

abundance Xc, nuclear energy generation rate per 
unit mass e n , neutrino emission rate per unit mass 

e v , and L r around the carbon-burning flame at se- -pig 3 The same as Fig 2 but for Xc = 02 

lected stages of the flame propagation. To display 
the variation of L r , the dotted lines are used when 
L r < and the long dashcd-lines when L r > 0. 
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Fig. 2. — Total nuclear energy generation rate L n 
(solid line) and neutrino energy loss rate L v (dotted 
line) are shown as functions of mf (the mass coor- 
dinate at the carbon burning flame) for Xc = 0.5. 
The phases with constant L D correspond to the phases 
during which the steady burning approximation was 
applied. Large fluctuations in L n at the ends of such 
phases are artifacts caused by inserting many grid 
points around the flame. 



Fig. 4. — Variation of the temperature distribution 
as the carbon-burning flame propagates inward. 
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Fig. 5. — Variation of the density distribution as the 
carbon-burning flame propagates inward. 



Fig. 7. — Evolution of temperature distribution be- 
fore (upper panel) and after (lower panel) the re- 
activation of flame. 




Fig. 6. — Changes in the position of flame (dashed Mr 
line) and thermal front (solid line) for Xc = 0.5. The 

dotted line indicates the thermal front for X c = 0.2, Fi §- 8 -~ Evolution of density distribution before (up- 

for which the flame is always very close to the thermal P er P anel ) and aftcr ( lower P anel ) thc re-activation of 

front. flame - 
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Fig. 9. — Propagation speed \dm{/dt\ (M Q y _1 ) ver- 
sus position mi of the flame for Xq = 0.2. Filled 
circles represent propagation speed evaluated from 
the propagation of the flame in evolutionary mod- 
els, while open circles represent the values obtained 
by interpolating Timmes et al.'s (1994) table for the 
density and temperature at the flame. 
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Fig. 11. — The time variations of L n and L v during 
the evolution phase shown in Fig. 10. 
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Fig. 10. — The upper panel shows L n and L v as func- 
tions of the mass coordinate of the flame for Xq — 0.2. 
Time runs from right to left. The lower panel shows 
temperature distribution near flames at various evo- 
lutionary stages. 



Fig. 12. — The time variation of the integrated nu- 
clear energy generation rate during the process of 
inward propagation of the helium burning shell for 
M = lO^Moy" 1 . 
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Fig. 13. — The mass coordinate at the helium burning 
shell versus time after the ignition. 
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Fig. 15. — Runs of temperature (upper panel) and the 
rate of gravitational energy release (e g , lower panel) 
around the temperature maximum at some selected 
stages during one cycle of a shell flash (from a peak 
of shell burning to the next one). 



Fig. 14. — The HR diagram showing evolutionary 
paths of accreting white dwarfs with M = 1.0 x 
lO" 6 M y- 1 (solid line) and M = 1.0 x lO^Afoy" 1 
(dotted line). (The accretion is terminated when the 
burning shell reaches the center.) 



14 



